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The Fermi surface in the pseudogap [1 —5] state of cuprates is highly unusual 
because it appears to consist of disconnected segments called arcs [6]. Their very 
existence challenges the traditional concept of a Fermi surface as closed contours 
of gapless excitations in momentum space. The length of the arcs in the pseu- 
dogap state was thought to linearly increase with temperature, pointing to the 
existence of a nodal liquid state below T*[7, 8]. These results were interpreted as 
an interplay of a d-wave pairing gap[9] and strong scattering[10]. Understand- 
ing the properties of the arcs is a pre-requisite to understanding the origin of 
the pseudogap[8, 11-20] and the physics of the cuprates. Here we use a novel 
approach to detect the energy gaps based on the temperature dependence of 
the density of states. With a significantly improved sensitivity, we demonstrate 
that the arcs form rapidly at T* and their length remains surprisingly constant 
over an extended temperature range between T* and T arc , consistent with the 
presence of an ordered state below T*. These arcs span fixed points in the mo- 
mentum space defining a set of wave vectors, which are the fingerprints of the 
ordered state that causes the pseudogap. 

Traditional methods for determining the presence of an energy gap in ARPES spectra 
[4-7] rely on a lineshape analysis such as a shift in the leading edge or a dip in symmetrized 
spectra[6]. In the cuprates however, these features are poorly defined (especially above 
T c ) because the spectral peaks are very broad in optimally and underdoped samples. This 
makes the detection of very small energy gaps and partial gaps (such as the pseudogap, 
which affects only part of the spectral weight) difficult or impossible to do. Even gapped 
ARPES spectra may appear to have a single peak at Ep after symmetrization if they are 
sufficiently broad. One of the most sensitive ways to detect of the opening of an energy 
gap is via the density of states (DOS(E'f)) at the Fermi energy. In a gapless state, this 
quantity is independent of temperature, because unlike the spectral function it does not 
depend on the electron lifetime. The opening of an energy gap at the Fermi energy leads 
to a decrease of DOS (.Bp)- Studying this quantity as a function of temperature therefore 
provides a very sensitive measure of a gap opening, which is quantitative and objective. To 
obtain momentum dependent information we use the area of momentum distribution curve 
(MDC) at the Fermi energy (Dmdc(-Ef)) which represents a contribution to the DOS(E'f) 
from a small slice of the Brillouin zone. This new approach requires very high quality 
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data and temperature stability, but at the same time it allows us to obtain momentum 
resolved information about the opening of an energy gap with significantly higher accuracy 
than previously possible. Details of the sample preparation and ARPES measurements are 
provided in the Supplementary Information. 

In Fig. 1 we demonstrate the above procedure using data at the node (gapless) and away 
from the node (gapped) from a Bi2201 sample and compare it to the traditional method 
of symmetrization. Fig. la shows the MDCs along the nodal cut (see inset of Fig. lb) 
measured at various temperatures from below the superconducting transition temperature 
(T c ) to above the T*. The MDC peak, which is observed at k F , broadens with increased 
temperature, as the electron lifetime shortens due to scattering. The area of the MDCs, 
Dmdc(-Ef), remains constant with temperature as expected for a partial contribution to the 
DOS in the absence of an energy gap (Fig. lb). This behavior changes dramatically away 
from the node, where an energy gap opens. Fig. lc shows the MDCs along a momentum 
cut close to the antinode (see inset of Fig. Id). At high temperatures above T*, the peak 
width increases, while the area remains constant. Below T*, the area of the peaks starts 
to decrease, indicating the opening of an energy gap. This is even more evident, when 
examining the Dmdc(-^f) as a function of temperature shown in Fig. Id. This quantity is 
constant at high temperatures, then decreases below T*. The reduction of D MDC (_E I F ) below 
T* signifies the opening of an energy gap - in this case the pseudogap (indicated by a blue 
arrow). We note that there are number of features in D MDC below the temperature at which 
the gap opens. They very likely contain important information related to the paring [21] 
and the superffuid density. They are, however, beyond the scope of this work and will be 
analyzed carefully in the future. 

For a comparison, we plot the symmetrized energy distribution curves (EDCs) at kp for 
several temperatures in Fig.le. This is the traditional way to detect the energy gap [6, 7]. The 
single peak present at high temperatures develops a dip on cooling, which was previously 
interpreted as a signature of a gap opening. It is clear that the temperature (T pea k) is 
considerably lower than the T*, because smaller gaps even if present do not always produce 
a dip in the symmetrized spectra. The traditional approach fails to detect small energy 
gaps and lacks the sensitivity necessary to reveal the real nature of the Fermi arcs. Careful 
investigation of Dmdc(Ef,T) is essential to achieve this goal. 

In order to study the temperature evolution of the Fermi arc, we carefully measured 
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Dmdc(Ef,T) for a number of Fermi momentum points. The sample orientation was ad- 
justed in each case, so the cuts were normal to the Fermi surface (see insets of Fig. 2). 
The left and right panels in Fig. 2 show the data obtained from optimally doped Bi2201 
(T C =32K) and Bi2212 (T C =93K), respectively (raw MDC data is included in the Supplemen- 
tary Information). The data in the top panels was obtained at the antinode, while the lower 
panels show data towards the node. The D M y) C (E f ,T) starts to decrease upon cooling at T* 
in the top four panels for 0° < <fi < 19°, indicated by blue arrow. For the remaining Fermi 
momentum points closer to the node, the Dudc(Ef,T) remains constant down to a much 
lower temperature. This signifies the simultaneous opening of the pseudogap for a range of 
Fermi momentum points up to < 19° at T*. The Fermi surface therefore collapses very 
abruptly into arcs just below T*. The lack of a variation in Dubc{Ef,T) for > 19° down 
to T arc demonstrates that the length of the arcs remains constant in this temperature range. 
On further cooling below T arc , D M dc(E f ,T) starts to simultaneously decrease over the re- 
maining segments of the FS. This signifies the second collapse of the remaining segments of 
the Fermi surface (arcs) below T arc into a very small arc or most likely the point node of 
a rf-wave paired state. We summarize these results in Fig. 3 by plotting the gap opening 
temperatures as a function of the Fermi surface angle and the Fermi arc length as a function 
of temperature, based on the data in Fig. 2. 

Figure 4 shows a schematic diagram demonstrating the process of forming the arcs and 
their subsequent collapse. Our most important finding is the double collapse of the Fermi 
surface on cooling first at T* and then at T arc . This result clearly contradicts the nodal liquid 
scenario [7, 8], where the Fermi arcs expand linearly with temperature. We emphasize that 
the traditional lineshape analysis underestimates the temperature at which the pseudogap 
opens. To confirm this, we repeated a traditional symmetrization analysis for our data from 
Fig. 2, and demonstrated an apparent linear expansion of the arcs with temperature (see 
Supplementary Information). The remaining open issue is whether a tiny Fermi arc or a 
single point exist very close to the node below T pair . At this time, due to the limits of the 
data accuracy we can only impose limits on the length of this arc to about 20% of the normal 
state Fermi surface. 

The Fermi momentum points at the tip of the pseudogap arcs coincide with the ends of 
the parallel segments of the FS rather than at the antiferromagnetic zone boundary (AFZB). 
The vectors connecting these points such as (0, 0.4), (1.4,0), (1.4, 0.4) are likely the key 
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to determining the nature of the ordered state underlying the pseudogap phenomenon in 
cuprates. Earlier STM studies [22, 23] proposed a scenario for the pseudogap based on a 
charge density wave with a (0, 0.4) ordering vector. The strong energy dependence of the 
checkerboard pattern [24] in STM is however inconsistent with the signatures of a classical 
charge density wave. The other two vectors however may be related to a nematic state 
reported by separate STM study [25]. So far we have not been able to identify a theory 
consistent with the above set of vectors, however we hope that this information will lead to 
the development of a correct theory of the pseudogap state. 
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FIG. 1: Determination of the partial density of state at Ep, Dmdc(£f) 5 from the ARPES data 
of optimally doped Bi2201 (T C =32K). a, MDCs measured along the nodal direction over a wide 
range of temperatures. Inset shows a dispersion image at the lowest temperature (12K). b, The 
temperature dependence of Dmdc(-Ef) estimated from the area of MDCs in a. The observed 
momentum cut and estimated area of MDC are demonstrated in the inset, c, Same data as in 
a, but measured along a momentum cut slightly off the antinode. Inset shows a dispersion image 
at the lowest temperature (12K). d, The temperature dependence of Dmdc(-Ef) estimated from 
the area of MDCs in c. The observed momentum cut and estimated area of MDC are shown in 
the inset. The pseudogap temperature (T*) is defined as the temperature where the ~Dmdc(Ef,T) 
starts to decrease on cooling, e, Symmetrized EDCs at k-p for several temperatures, obtained from 
the data of c. The T pca k is defined as the temperature where two peaks in the spectrum merge to 
one peak at elevated temperatures. 
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FIG. 2. Angle dependence of the partial density of state, D mdc (-Ef)- a— i, for optimally 
doped Bi2201 with T C =32K. j-r for optimally doped Bi2212 with T C =93K. The top panels 
correspond to the antinodal direction. Lower panels show data along cuts approaching the 
nodal direction. The inset of each panel shows the FS and the location of the cut, normal to 
the FS. The pseudogap temperature (T*) and the temperature above which the arc exists 
(T arc ) are indicated by blue and red arrows, respectively. 
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FIG. 3: Momentum dependence of the temperatures at which gap opens and the length of the 
Fermi arc determined from ARPES data shown in Fig. 2. a, Gap opening temperatures as a 
function of Fermi surface angle for optimally doped Bi2201. b, Temperature dependence of the 
arc length. Abrupt changes in the length of the arcs occur at T arc and T* . Their length remains 
constant for T arc <T <T*. c— d, same as a and b, but for optimally doped Bi2212. 
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FIG. 4: Schematic diagram of the temperature evolution of the Fermi surface in cuprates revealed 
by the new method, a— c, The Fermi surface for three key temperature ranges. The antiferro- 
magnetic zone boundary (AFZB) is marked with a black dashed line. The black dotted lines are 
the guide to eyes for parallel segments of Fermi surface, d, The Fermi surface for the Bi2212 and 
Bi2201 determined from our ARPES data. The key q vectors connecting the nested sections of the 
FS are indicated with arrows. The black dotted lines a the guide to the eye for parallel segments 
of Fermi surface. 
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I. SAMPLES AND EXPERIMENTAL METHOD OF ARPES 

Optimally doped Bi 2 Sr 2 CaCu20 8+( 5 (Bi2212) single crystals with T C =93K (OP93K) and 
(Bi,Pb) 2 (Sr,La) 2 Cu0 6+5 (Bi2201) single crystals with T C =32K (OP32K) were grown by the 
conventional floating-zone (FZ) technique. We partially substituted Pb for Bi in Bi2201 to 
suppress the modulation in the BiO plane, and avoid contamination of the ARPES signal 
with diffraction images due to the superlattice: the outgoing photoelectrons are diffracted 
at the modulated BiO layer, creating multiple images of the bands and Fermi surface that 
are shifted in momentum. The modulation-free samples enable us to precisely analyze the 
ARPES spectra. 

ARPES data was acquired using a laboratory-based system consisting of a Scienta 
SES2002 electron analyzer and GammaData Helium UV lamp. All data were acquired 
using the Hel line with a photon energy of 21.2 eV. The angular resolution was 0.13° and 
~ 0.5° along and perpendicular to the direction of the analyzer slits, respectively. The 
energy corresponding to the chemical potential was determined from the Fermi edge of a 
polycrystalline Au reference in electrical contact with the sample. The energy resolution 
was set at ~ lOmeV - confirmed by measuring the energy width between 90% and 10 % 
of the Fermi edge from the same Au reference. Custom designed refocusing optics enabled 
us to accumulate high statistics spectra in a short period of time with no sample surface 
aging from the absorption or loss of oxygen. Special care was taken to purify the helium gas 
supply for the UV source to remove even the smallest trace of contaminants that could con- 
tribute to surface contamination. Typically no changes in the spectral lineshape of samples 
were observed in consecutive measurements performed over several days. We constructed a 
sample manipulator with the tilt and azimuth motions mounted on a two stage closed cycle 
He refrigerator. To measure the partial density of states along a direction perpendicular 
to the FS we controlled the sample orientation in-situ. Measurements were performed on 
several samples and we confirmed that all yielded consistent results. 

II. NORMALIZATION OF ARPES SPECTRA 

To reduce the noise level when collecting data over an extended temperature range, we 
normalized the band dispersion map (see Fig. SI e) at each temperature, to the whole 
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area of the symmetrized energy distribution curve (EDC) at the kp point in the range of 
-0.4eV < E < 0.4eV. Here we demonstrate that a particular normalization scheme does not 
affect the results. Fig. Sf a and c show the symmetrized EDCs at the k F point near the 
antinode measured at various temperatures. The spectra in panel a are normalized to the 
whole area of each spectrum in the range of -0.4eV to 0.4eV, while spectra in panel c are 
normalized to the intensity within ±50 meV centered at -0.35eV. Please note that in the 
data normalized over the wide energy range (panel a), the intensity at higher energies (e. g. 
-0.35eV) is almost constant for all temperatures. It is, therefore, expected that the results 
obtained by the two different normalization schemes are consistent. The data normalized 
using these two procedures at each temperature was used to extract MDC area at Ep (i. e. 
Dmdc(-Ef) vs T) and is shown in Fig. IS f. As expected, the temperature dependence of 
this quantity is very similar for the two normalization schemes. The two curves are almost 
identical with the same onset of the reduction in Dmdc(-^f) on cooling at T*=180K (the 
pseudogap temperature). This allows us to conclude that the results we present, including 
the determination of the gap opening temperatures, do not dependent on the choice of a 
normalization scheme. 

The intensity of spectral function is expected to be zero far from the location of the 
band. However, the ARPES intensity at E F in such a momentum region is usually finite. 
This background signal is due to photoelectrons that have lost their momentum information 
during the photoemission process via elastic scattering[l]. We subtract this intensity from 
the experimental momentum distribution curves (MDCs) before estimating the MDC area 
for each temperature. We have checked the validity of this treatment by performing the 
same analysis for many data sets from different samples. The intensity due to the elastic 
scattering varies slightly from sample to sample because elastic scattering is very sensitive 
to the condition of the sample surface. The behavior of D M dc(E f ) (area of MDC) estimated 
before and after such a procedure is nevertheless consistent for all samples. The checks 
described above validate all steps of our analysis of the ARPES data. 
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III. COMPARISON OF THE RESULTS FROM QUANTITATIVE ANALYSIS 
WITH LINESHAPE ANALYSIS 



In Fig. S2 and S3 we use our high-quality data to demonstrate the inadequacy of the 
traditional lineshape analysis whereby the spectral peak in the symmetrized data is tracked. 
This method results in an underestimation of the temperature at which the gap closes. 
Figure S2 and S3 show the symmetrized energy distribution curves (EDCs) measured at 
various temperatures from Bi2201 and Bi2212 samples, respectively. From left to right, 
we plot the data from the antinode to the node. At elevated temperatures, the two peaks 
in the spectra merge into one (a red curve) at a characteristic temperature (T peak ). It is 
easily seen that the T pea k monotonically decreases toward the node. As a result, the portion 
of the Fermi surface where the spectra have a peak at Ep emerges centered at the node, 
and it expands toward the antinode with increasing temperature. In the inset of Fig. S2 
and S3 we estimate the length of these "artificial" Fermi arc as a function of temperature. 
Such plotted arc lengths are similar to an earlier report [2]. The length of arcs increases 
almost linearly with temperature above T c up to T*. We emphasize that at momentum 
points, where the traditional lineshape method determines the gap to be closed, the partial 
density of states is still increasing with temperature. This signifies that the energy gap is 
still present, but is not detected using symmetrization and is most likely due to thermal 
broadening or a rapid increase of scattering above T c . The biggest advantage of utilizing 
the partial density of states at E F , D M dc(E f ,T), is that it is not affected by the thermal 
broadening and is very sensitive to presence of even a very small energy gap. From the 
reduction in the D M p, c (E F ,T), one can estimate the gap opening temperature much more 
precisely than with the traditional approach. 

IV. RAW DATA USED IN ANALYSIS SHOWN IN FIG. 2. 

In Fig. S4 a-i and S5 a-i, we present the MDCs along many directions of Bi2201 OP32K 
and Bi2212 OP93K, respectively, used to estimate the partial density of states, D M dc(-Ef), 
(shown in Fig. 2 of the main text). In the inset of each panel we plot the band dispersion 
image for the lowest temperature (12K). The panels j-r of Fig. S4 and S5 show the same 
data as in Fig. 2. In each panel, the gap closing temperature estimated from the spectral 
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peak position (T pcak ) and the density of state (T D qs) are indicated. 
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Fig. SI. Verification of normalization procedure, a, Symmetrized EDCs at the kp close to 
the antinode [see inset of f] measured over a wide temperature range from below T c to above 
T*. All spectra are normalized to the total area integrated from -0.4eV to 0.4eV. b, MDCs 
at various temperatures along the momentum cut indicated in the inset of f, extracted from 
the band dispersion maps normalized to the total area of symmetrized EDCs at kp for each 
temperature, c, Same spectra as in a normalized to the intensity within a narrow energy 
window (±50meV) centered at -0.35eV. d, MDCs at various temperatures, extracted from 
the dispersion map normalized to the intensity of spectrum at kp within a narrow energy 
window (±50meV) centered at -0.35eV. e, Band dispersion map at the lowest temperature 
(12K). f, Two sets of D MD c(-Ef) v s T curves [area of MDC], obtained using two different 
normalization schemes. The inset shows the estimated area of MDC and the Fermi surface 
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with the momentum cut measured. 
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Fig. S2. The "artificial" Fermi arc of Bi2201 OP 32K obtained by tracking spectral peak. 
In each panel, the symmetrized EDCs measured over a wide temperature range from below 
T c to above T* are plotted. From the left (a) to the right (i) panel, the measured hp point 
changes from the antinode to the node. The Fermi angle of the measured hp point [0: see 
inset of j] is described on the top of each panel. The spectrum where the two peaks merge 
to one peak at elevated temperatures is indicated with a red color, j, The length of Fermi 
surface ("artificial" Fermi arc) where one peak structure is observed is plotted as a function 
of temperature. 
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Fig. S3. The "artificial" Fermi arc of Bi2212 OP 93K obtained by tracking spectral peak. 
In each panel, the symmetrized EDCs measured over a wide temperature range from below 
T c to above T* are plotted. From the left (a) to the right (i) panel, the measured kp point 
changes from the antinode to the node. The Fermi angle of the measured kp point [0: see 
inset of j] is described on the top of each panel. The spectrum where the two peaks merge 
to one peak at elevated temperatures is indicated with a red color, j, The length of Fermi 
surface ("artificial" Fermi arc) where one peak structure is observed is plotted as a function 
of temperature. 
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Fig. S4. Direction dependence of the partial density of state at E F , D MDC (E F ), for the 
optimally doped Bi2201 with T C =32K. a-i, MDCs at various temperatures from below T c 
to above T*. The top panel a corresponds to the antinodal cut. Toward the bottom i, the 
momentum cut measured is shifted toward the node. Each inset shows the band dispersion 
map, from which the MDC is extracted, at the lowest temperature (12K). j— r, The tempera- 
ture dependence of the partial density of state, D MDC (i? F ), estimated from the area of MDC 
in a-i. The inset of each panel shows FS with a measured momentum cut, that is normal to 
FS. The gap opening temperature (T D os) estimated from the reduction of Dmdc(-Ef) and 
the temperature at which the symmetrized EDC has a peak at E F are indicated with a blue 
and red dotted line, respectively. 
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Fig. S5. Direction dependence of the partial density of state at E F , D MDC (_E F ), for the 
optimally doped Bi2212 with T C =93K. a-i, MDCs at various temperatures from below T c 
to above T*. The top panel a corresponds to the antinodal cut. Toward the bottom i, the 
momentum cut measured is shifted toward the node. Each inset shows the band dispersion 
map, from which the MDC is extracted, at the lowest temperature (12K). j— r, The tempera- 
ture dependence of the partial density of state, D MDC (i? F ), estimated from the area of MDC 
in a-i. The inset of each panel shows the FS with the measured momentum cut, that is nor- 
mal to FS. The gap opening temperature (T D os) estimated from a reduction of Dmdc(-Ef) 
and the temperature at which the symmetrized EDC has a peak at E F are indicated with a 
blue and red dotted line, respectively. 
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